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Abstract 

We analyze the tensor perturbations of flat thick domain wall branes in f(R) gravity. Our results indicate that 
under the transverse and traceless gauge, the metric perturbations decouple from the perturbation of the scalar field. 
Besides, the perturbed equation reduces to the familiar Klein-Gordon equation for massless spin-2 particles only when 
the bulk curvature is a constant or when f(R) = R. As an application of our results, we consider the possibility of 
localizing gravity on some fiat thick branes. The stability of these brane solutions is also shortly discussed. 
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1. Introduction 

The trapping of various matter fields on four-dimensional domain walls has been successfully realized by using 
both non-gravitational or gravitational methods HI 01- Recently, domain walls in five-dimensional space-time (called 

J~~| |~~| 

the bulk) attract renew attentions from the physical community after Randall and Sundrum yl |4J pointed out that 
massless four-dimensional graviton can be realized on a thin wall (called the brane) if the extra dimension is large and 
warped properly. Randall and Sundrum assume that all the observable matter fields are confined on the brane while 
gravitons can propagate in the bulk. For this reason, their model is also known as the (thin) braneworld scenario, or 
the RS model for short. The braneworld scenario is very important, because it, to some extent, solves some of the long 
existed problems such as the hierarchy problem, and the cosmological constant problem, etc., for reviews see 

However, the thin braneworld model suffers a drawback: the whole theory is singular at the location of the brane, 
for example, the bulk curvature diverges at there and a junction condition must be introduced. Besides, in the original 
braneworld scenario we are not clear how does our world (as a brane) formed. Such problems are solved in the so 
called thick braneworld models, in which gravity couples with a background scalar field j7l 4l ill . The formation of the 
domain wall is ascribed to the non-linear property of the gravitational system. 

In general relativity, four-dimensional gravity has been successfully realized in some thick braneworld models 10- 
fioll . The trapping of various kinds of matter fields on single or multi branes are also discussed for both thin and thick 
branes I12h31I1 . 

However, for some reasons, we have to take the contributions from the higher order curvature terms into consid- 
eration. On one hand, the adding of higher order curvature terms might render general relativity renormalizable [32]. 
On the other hand, higher-order curvature invariants also appear in the effective low-energy gravitational action of 
more fundamental theories, such as the string theory 13311 . To prevent the theory from the spin-2 ghosts, and the well- 
known Ostrogradski instability 13411 . the higher order curvature terms are usually introduced as the Gauss-Bonnet term 
or an arbitrary function of the curvature, namely, f(R). Both of these two modified gravity theories were applied in 
discussing a wide range of issues in cosmology and higher energy physics, for details, see B351437I1 and references 
therein. 

It was shown that the introduce of the Gauss-Bonnet term usually imposes no impact on the localization of gravity 
on the brane However, in f(R) gravity things are more complex, since the dynamical equations are of 
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fourth-order. Based on the fact that f(R) gravity is conformally equivalent to a second-order gravity theory [48], 
some thin braneworld models have been constructed in the lower order frame | U 52 ] by introducing a conformal 
transformation. We hope this method still valid for the discussions of thick branes. However, as stated in 115 311 . the 
method used in refs. 149145211 would lead to ambiguous when a background scalar field is introduced. 



For this reason, the thick brane solutions usually were found directly in the higher order frame 154 15511 . In [54J], 
with a background scalar field, the authors offered us some analytical thick brane solutions in both constant curvature 
spaces and more general space-time. While in 15511 . the authors numerically discussed some thick brane solutions 
with pure gravity. In fact, the authors of lf55ll identified the contribution of higher order curvature term f(R) — —aR" 
(a > and 1 < n < 2) with an effective "matter" source. The solutions were obtained by analyzing the existence of 
the fix points. Besides, the trapping of complex scalar field on the brane solutions was also discussed in 15511 . 

An ideal thick brane solution should be smooth, stable and possible to localize gravity and various kinds of matter 
fields. Whether the four-dimensional gravity can be reproduced on the brane solutions in 15415511 is still unclear. In 
order to address the issue of localizing gravity on thick /(/?)-brane, we analyze the tensor perturbations for a particular 
model. Finally, we apply our results on some of the solutions found in f54ll . 

We organize this letter as the follows: In the next section, we set up our model and give the dynamical equations. 
In section[3]we discuss the tensor perturbations of the model. The localization of four-dimensional gravity on some 
thick branes is investigated in section[4] We focus mainly on the solutions for which the bulk curvatures are constants. 
For such solutions the perturbed equations reduce to the Klein-Gordon equation for massless spin-2 particles. Our 
summary is given in section|5] 



2. The model 

We start with the action 



S = J ctxdysf—g 



JL /CR) _ l -d M 4>d M $ - V{<f>) 



(1) 



where <p is a background scalar field which generates the brane. V(<f>) describes a self-interacting potential for the 
scalar field. The gravitational coupling constant k 2 = 87rG5 with G5 the five-dimensional Newtonian constant. In- 
dices M,N ■ ■ ■ = 0, 1 , 2, 3, 5 and fi, v ■ ■ - 0, 1 , 2, 3 are always applied to denote the bulk and the brane coordinates, 
respectively. 

For simplicity, we consider the static flat braneworld scenarios with the metric 

ds 2 = e 2A(y) r lMV dx^dx v + dy 2 , (2) 

and the scalar field is assumed to be a function of the extra dimension y = x 4 , i.e., <p - <P(y)- For system (Q]i-©, 
Einstein equations are 

f(R) + 2f R (4A' 2 + A") - 6fgA' - 2% = k\{$' 2 + 2V), (3) 



and 



8/« (A" + A' 2 ) + 8JJA' - f(R) = k\(<S>' 2 - 2 V), (4) 



where the primes represent derivatives with respect to the coordinate y and fg = df(R)/dR. The equation of motion 
for the scalar field is given by 

dV 

4A>' + </>" - — = 0. (5) 

0(p 

Obviously, these equations contain higher-order derivatives of the coordinates. Usually, it is very hard to solve these 
equations analytically, needless to say to analyze the feedbacks of these equations for the perturbations from both 
the metric and the scalar field, because, in general, the perturbed equations are coupling equations with higher-order 
derivative terms. 
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3. Tensor perturbations 

We consider the following metric perturbations: 

ds 1 = e 2A(y \r] MV + h llv )dx^dx v + dy 2 , (6) 



or, in another form 



with 



gMN - gMN + k-gMN, (7) 



8mn = \ q l I , &g MN = I o I ( 



the background metric and the metric perturbations, respectively. Here = h MV (x p , y) depend on all the coordinates 
Obviously, Ag$M = 0, which means we consider only tensor perturbations. According to the relation g NP gpM = 6$- 
one obtains the inverse of AgMN, i-e., Ag MN . We keep only the first order term, i.e., Ag MN( - r> and denote it as 



where W = rf A T] vp h^ p is raised by rf v . We always use 6X to denote the first order contribution of the perturbations 
to an arbitrary quantity X. The perturbation of the scalar field is assumed to be of first order and is denoted by 
6<p = y). 

Denoting a(y) = e A(y \ we immediately obtain the following relations: 

5R I1V = -i(D (4) /i^ v + d M d v h- d v do-h* - d^da-h^) -laa'h'^ 



3h uv a' 2 - ah uv a" - 



„ _ a Ky _ aVuva'ti 

/(!'" II//,,, It - — , 



6R . w = _« + ^_^-,,». ,10, 

a 1 a 1 a 

Here D (4) = Tf v d M d v , is the four-dimensional d' Alembert operator, and h = rf v h IJV is the trace of the tensor perturba- 
tions. 

Obviously, if 

h = = d li h" v , (11) 

only SR^ remains non-zero. The condition (TTTT i is called the transverse-traceless gauge which can largely simplify 
the perturbed equations. 

In f(R) gravity, the Einstein equations are 

Rmn/r - -jgMNfiR) + (8mnO (5> - V M V N )f R = k]T mn , (12) 

with D (5) = g MN V m^n the five-dimensional d' Alembert operator. As the perturbations are considered, the feedback 
from Einstein equations reads 

6Rmn/r + RmnIrrSR - 2^8MN.f(R) - ^Smn/rSR 
+ 6(gMNa (5) f R )-6(V M V N f R ) = K 2 5 6T MN . (13) 
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We need only to calculate the second line of the above equation. Note that 

gMNa (5) f R = gMNg AB (^A^BfR), 



we have 



£(VmV w /r) = (d M d N - TM N d P )(f RR SR) - 5T P MN d P f R , 

6{,g MN U {5) f R ) = SgMN^ 5) fR+gMN5g AB {V A V B f R ) 
+ gMNg AB 6(V A V B f R ). 



The fluctuations of the energy-momentum 

Tmn 

are given by 



V M <pV N 4> - -gMNg 1 V A <^V fi - gMNV((p) 



ST, 



U ~ -fit. „2„ J in . 2 9V~. 



fiv 



ST 5fI 



<p' h^, - a ri^'P - a Vh^y - a rj^- 
.,. dV~. 
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Note that under the transverse and traceless gauge, 6R = 0, and eqs. ( fT6b reduce to 

S(V M V N f R ) = -6 P M 6%6T 5 pcr f R , 
6(g MN D^f R ) = 5 p M 5«5g ptT n {5) f R . 



Therefore, under this gauge, we have 



6(g M ND (5) f R ) - S(V M V N f R ) 



= *» 2 



The perturbed Einstein equations ( fTJb reduce to 



SRmn/r ~ -^SgMNfiR) 



By plugging ( TTOt and ( TT~8T > into (fJTJ, we obtain the (ft, v) components of the perturbed Einstein equations 



1 a 2 h" v ^ 
._□( )/j^ v _ Ih^a' 1 - laah'^ - ah^a" 



.h 



1 



-a h^yfiR) + a 



V V'—fk +/r 



1 



^ l livf R 



2 ?liJY<t>'<i>' - c^Vh^y - a 2 77 /jv ^0 



Note that the (/i, //) components of Einstein equations ( fT2l > is 



f(R) + 2f R 



which is exactly Eq. (0. We can simplify Eq. ( l22l as 



1 a 2 h" x 



= -cfijftytp'ip' - a 2 ^ 



(24) 



Contracting the above equation with rf Y one can proof that + = 0. Therefore, the (/j., v) components of the 
perturbed Einstein equations read 



or, equivalently, 



With the coordinate transformation 



- 2 D (4 V + 4-h' + h")f R + h'f R = 



□ (5) ^V = T-dyh^y. 

JR 



dz — a l dy, 



we can rewrite the perturbed equation ( t26l > as 



a /« 



V = o. 



(25) 



(26) 



(27) 



(28) 



Consider the decomposition h^ v (x p ,z) = (a 3/2 f R ll2 )e llv (x p )>l/(z), and ask e /iv (jc p ) satisfies the transverse and traceless 
condition rf v £ flv = = S^e/, we would have a Schrodinger equation for if/(z): 



[d 2 - W(z)] iff(z) = -m 2 iff(z), 



with the potential W(z) given by 



W{Z> A a 2 2 a 2 af R 4 f 2 2 f R 



To be more explicit, one can also factorize the Schrodinger equation d29l as 

i^iz) = -m 2 if/(z), 



2 a 2 f R 



() ■ 3 d : a | 1 d,/ R 
2 a 2 



(29) 



(30) 



(31) 



which indicates that there is no gravitational mode with m 2 < 0. Therefore any solution of the system <JTJ-([2} is stable 
under the tensor perturbations. The zero mode (if exists) takes the form 



^ \z) = N a^(z)fr(z) 



3/2/ 



d/2/ 



(32) 



with No the normalization constant. 

These results indicate that as the transverse and traceless gauge is taken, the perturbation of the scalar field de- 
couples from the metric perturbations. In the case of general relativity, Eq. (l26l > reduce to the five-dimensional Klein- 
Gorden equation for the massless spin-2 gravitons. However, for an arbitrary form of f(R) and non-constant curvature 
R, the equation for is largely different from the massless Klein-Gorden equation. Fortunately, the perturbed equa- 
tion always remains second order due to the introducing of the transverse and traceless gauge. Now let us see a simple 
application of our results on some solutions given previously in 1 54] . 
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4. Applications 



In 115411 . the authors gave us brane solutions for both constant and variant curvature cases. However, the solution 
for later case contains a singular point, and therefore is not regularized. So we would like to consider only the constant 
curvature case. The corresponding Einstein equations in this case reduce to some second ones: 



f(R) + 2f R 



b' 2 + 2V), 



-8/«— -f(R) = 4(<f>' I -2V). 
a 



Constraining R as a constant, the solution for warp factor is uniquely determined [5 

for R = 0, 



a(y) 



cos 



for R = 20k 2 > 0, 



cosh 2/5 ^kyj , for R = -20k 2 < 0. 



(33a) 
(33b) 

(34a) 
(34b) 

(34c) 



Note that the warp factor given in ( 134at is not smooth, it contains a cusp at y = 0. This cusp of the warp factor leads 
to at most a (5-function in the second-order Einstein gravity, which can be explained as the appearance of a thin brane. 
However, in the fourth-order f(R) gravity, such cusp will introduce derivatives of 6(y), such as 6'(y) and S"(y), and 
cross terms of them. It is still a problem of how to deal with such terms in brane theory. For this reason, we consider 
only solutions for dS 5 space where R > 0; and for AdS 5 space where R is a negative constant. 

One can easily proof that for dSs and AdS 5 spaces, eqs. (l33l support the following non-trivial solutions: 



For dS 5 and fn>0 




V(cf>) 



+ _ /^^-arctanh(sin(5Ary/2)), 
5/c 



V.-^sintf 

4k 5 




(35) 
(36) 



with V] 



2f(R)-25k 2 f R 



We consider the interval —£r < y < within which the warp factor ( I34bl ) is regularized. 



For AdS 5, and f R < 



4> 

V(4>) 



6|/j?l 



V 2 + 



A/ 5* 

9k 2 \f R 



arctan (sinh(5Ary/2)) , 



4k 5 



sin 



5k 2 5 



\6|/*| 



(37) 
(38) 



with v 2 = m - 2 4m. 

For the warp factors given in (I34bb and (I34cl ). the coordinate transformation d27l i cannot be integrated out ana- 
lytically. However, we can find the numerical relation between z and y, see figure Q] Note that for the case of dS 5 
space, 



z(y) 



cos 5 i-kyjdy = — - 

Jo V 2 > 5fer (4) 



(39) 
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Figure 1: Figure of z(y) for AdS 5 space (left panel) and for dS$ space (right panel) with k = 2/5. 



where T(fi) — f Q l e 'dt is the Euler gamma function. Therefore, z(y) is a bounded function of y. 
For constant curvature spaces f' R = 0, the perturbed Einstein equations reduce to 

□ (5) V = 0. (40) 

This is the familiar Klein-Gordon equation one obtains in general relativity |^,@]. To analyze the possibility of 
localizing four-dimensional gravity on the brane, we simply solve the shrodinger equation d29b with the following 
potential: 

3 a' 2 3 a" 

4 a z 2 a 

The corresponding zero mode takes the form 

iff (0) (z) oc a 3/2 (z). (42) 

As shown in figure [2] for the case of AdS*, space-time, the potential W(z) is positive everywhere and diverge at 
Z = ±00. Such potential supports only bounded and discrete KK states. There is no continuum spectrum. Thus, the 
solution itself is stable under tensor perturbations. However, the zero mode dose not exist, as a consequence, the 
four-dimensional massless graviton cannot be localized on the brane. 

The potential W(z) for dS$ space, as shown in figure [2j is a bounded function. Both the potential W(z) and the 

zero mode (see figre[3} vanish at the boundary z(y - +^) = ± l ^ . That means the solution in the case of dSs is 
stable; in addition, the four-dimensional graviton can be localized on the brane. 



5. Conclusions 

To sum up, we considered the fluctuations from both the metric and the scalar field around the flat thick f(R)- 
branes. It turns out that the perturbation from the scalar field decouples from the tensor part of the metric perturbations 
when the transverse and traceless gauge is considered. The propagation of the metric perturbations in the bulk is not 
described by the simple massless Klein-Gordon equation any more, except f(R) = R or R is a constant. As an 
application of our results, we studied the stability of some simple solutions given previously. These solutions were 
found by constraining the bulk curvatures as constants. Among these solutions, the one for R = is problematic 
because at the location of the brane the metric poses a cusp, which would lead to the problem of divergence in the 
fourth-order f{R) gravity. The analysis of the solution in case of AdS 5 indicates that the solution is stable, there are 
infinite discrete massive KK states. While for the case of dS 5, the solution is stable, and the normalizable zero mode 
does exist. Therefore, the four-dimensional massless graviton can be localized on dSs brane. 
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Figure 2 : Figure of W(z) for AdS s space (left panel) and for dS 5 space (right panel) with k = 2/5. 
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z 

Figure 3: The zero mode if/^Hz) for the case of dS; space with k = 2/5. 



The results we obtained in this letter are valid for any solution of the system QJ-©. In the present letter, we 
confined our discussions only on the simplest case, i.e., R =const.. However, the constant curvature spaces are rather 
special, because in this case the equations of motion of the tensor perturbations are independent of the form of /(/?). 
Thus it is natural to ask what would be different if the curvature is variant. To answer this question we have to firstly 
find a well behaved thick /(/?)-brane solution which is regularized, stable, and analytical (rather than numerical). 
We also hope that the localization of gravity and the trapping of bulk matters are guaranteed by these solutions. 
Unfortunately, there is, to our knowledge, no such solution has ever been reported. In one of our recent works J56)], 
we have found that at least in squared curvature gravity (where f(R) cc R 2 ), there exists a thick domain wall solution 
which possess nearly all the properties we are searching for. It seems that for the case of variant curvature, the thick 
/(/?)-brane solutions contain more interesting features. 
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